arXiv:1501.03585vl [cond-mat.supr-con] 15 Jan 2015 


Multiband superconductivity in Ta4Pd3Tei6 with anisotropic gap structure 

Wen-He Jiao,^’El Yi Liu,^ Yu-Ke Li,^ Xiao-Feng Xu,^ Jin-Ke Bao,^ 

Chun-Mu Feng,2 S. Y. Li,4.5 Zhu-An Xu,^’^ and Guang-Han Cao^’^.Q 

^School of Science, Zhejiang University of Science and Technology, Hangzhou 310023, China 
^Department of Physics, Zhejiang University, Hangzhou 310027, China 
^Department of Physics, Hangzhou Normal University, Hangzhou 310036, China 
^ State Key Laboratory of Surface Physics, Department of Physics, 
and Laboratory of Advanced Materials, Fudan University, Shanghai 200433, China 
^Collaborative Innovation Centre of Advanced Microstructures, Nanjing University, Nanjing 210093, China 

(Dated: January 16, 2015) 

We carried out the measurements of magnetoresistance, magnetic susceptibility and specific heat 
on crystals of the low-dimensional transition metal telluride Ta 4 Pd 3 Tei 6 . Our results indicate that 
Ta4Pd3Tei6 is an anisotropic type-II superconductor with the extracted Ginzburg-Landau parameter 
i^Gh = 84. The upper critical field Hc 2 {T) shows a linear dependence at low temperature and the 
anisotropy of Hc 2 {T) is strongly T-dependent, both of which indicate a multiband scenario. A 
detailed analysis reveals that the electronic specific heat Gei(r) can be consistently described by a 
two-gap (s-|-d waves) model from the base temperature T/Tc ~ 0.12 up to Tc. Our data suggests 
multiband superconductivity in Ta4Pd3Tei6 with anisotropic gap structure. 

PACS numbers: 74.70.-b, 74.78.-w, 74.25.Op, 74.25.Bt 


Exploring exotic pairing mechanisms in new super¬ 
conducting compounds have always been one of the 
most attractive issues in condensed matter physics. 
Although the Bardeen-Cooper-Schrieffer (BCS) theory^ 
provides excellent explanations for the conventional 
phonon-mediated s-wave superconductor, the succes¬ 
sive emergence of unconventional superconductivity, 
claimed to be in either the singlet s±-wave^ and d- 
wave pairing statesf^ii^ or the triplet p-wave^ and 
/-wave pairing statesji presents great challenges for 
physicists to uncover the variously exotic non-phonon- 
mediated pairing mechanisms. The spin or charge 
fluctuations, which always appear in low-dimensional 
structures, were often regarded as the pairing glues for a 
large number of unconventional superconductors, such 
as the cuprates^ iron-based superconductors^^ quasi- 
one-dimensional (QID) superconductorS)2ii£ and heavy- 
fermion superconductors !^ Among them, multiband 
superconductivity and symmetry-imposed nodes were 
commonly observed as evidenced by a variety of tech¬ 
niques including thermodynamic approaches In this 
regard, the gained information on the gap symmetry has 
significant meanings in understanding the novel pairing 
mechanisms. 

Recently, we discovered bulk superconductivity with 
Tc = 4.6 K in a transition metal telluride Ta 4 Pd 3 Tei 6 '^ 
This material has a layered structure with QID charac¬ 
teristics consisting of PdTe 2 chains, TaTea chains and 
Ta 2 Te 4 double chains along the crystallographic 6-axis. 
The (103) planes, as illustrated in the inset of Fig. 1(c), 
are formed by the connection of the chains mentioned 
above along the [301] direction. These structural features 
are well reflected in the morphology of the crystals 
as shown in the inset of Fig. 1(a). For simplicity, 
hereafter we define a*-axis as to be parallel to the 
[301] direction and c*-axis as to be perpendicular to 


the (103) plane. Compared with layered chalcogenide 
superconductors NbaPd^SeT^ and M 2 PdxQ 5 {M = 
Nb and Ta, Q = S and Se)i^“— with similar QID 
structures, Ta 4 Pd 3 Tei 6 is a stoichiometric compound 
with flat two-dimensional sheets and it shows stronger 
electron-electron interactions Thermal conductivity 
measurements demonstrate the presence of nodes in 
the superconducting gap by the evidences of a large 
residual kq/T in zero field and the rapid increase of 
kq/T in magnetic field, mimicking the behavior of d- 
wave cuprate superconductor Tl-2201.— Furthermore, 
a temperature-pressure superconducting dome was as 
well observed, leading to the expectation that certain 
density-wave orders may neighbour the superconducting 
dome. According to the result of electronic structure 
calculations made by Singhj2^ any nearby density wave 
instability, if exists, would be a charge-density-wave 
(CDW) but not a spin-density-wave (SDW) as the 
transition metal contribution to the density of states 
is rather small. Besides, an s-wave pairing state from 
phonons associated with the Te-Te p bonding and a Fermi 
surface associated Te p bands are argued to be the most 
possible case. To explain the results of the thermal 
conductivity, the author also proposes the likelyhood 
of Ta 4 Pd 3 Tei 6 being a clean multi-gap superconductor 
where the gap ratio is significantly largei^ However, 
more recently, the low-temperature scanning tunneling 
spectroscopy (STS) measurements show the evidence for 
the presences of anisotropic superconducting gap with 
gap minima or even nodej^ Thus, the issue of the 
gap symmetry in Ta 4 Pd 3 Tei 6 is worthy to be addressed 
further. 

In this Rapid Communication, we investigate the 
superconducting properties of the Ta 4 Pd 3 Tei 6 crys¬ 
tals via systematic magnetoresistivity (MR), magnetic 
susceptibility and specific heat measurements. The 
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FIG. 1. (Color online) Magnetoresistivity of the Ta 4 Pd 3 Tei 6 crystal for (a) H \\ a*, (b) H \\ b and (c) H \\ c*. (d) The extracted 
superconducting upper critical field Hc 2 {T) for different field orientations. Werthamer-Helfand-Hohenberg (WHH) fitting is 
made to the data for H || b, represented by the red dashed lines. The black dashed lines show linear fittings. The inset of (a) 
shows the image of the crystal under an optical microscope, from which the crystallographic directions can be easily identified. 
The crystal structure of Ta 4 Pd 3 Tei 6 viewed perpendicular to the (103) plane is plotted in the inset of (c). The inset of (d) 
shows the anisotropy of Hc 2 {T). 


results indicate that Ta 4 Pd 3 Tei 6 is an anisotropic type- 
II superconductor. The upper critical fields i/c2(0) were 
estimated to be 4.2, 9.6 and 3.3 T for fields applied along 
the three axes, a*, b and c*, respectively. We observed 
strong T-dependent anisotropy of Hc 2 {T) and the linear 
increase i7c2, both of which suggest the multiband 
superconductivity. The zero-field heat capacity data can 
be best described by a two-gap model with anisotropic 
gap structure. In addition, the electronic specific heat 
coefficient in the mixed state, 7 e(i?), exhibits a nonlinear 
behavior. Thus, our results suggest two energy gaps 
with anisotropic gap structure are associated with the 
superconductivity in Ta 4 Pd 3 Tei 6 . 

Single crystals of Ta 4 Pd 3 Tei 6 were grown by a self-flux 
technique as previously describedJ^ High quality of as- 
grown crystals was confirmed by X-ray diffraction. To 
insure the data obtained more reliable, we elaborately 
selected crystals for the measurements, all of which 
having shiny surfaces without any discernible solvent 
attached. MR measurements were carried out by a stand 
four-probe technique with current 1 = 2 mA applied 
along the 6-axis. The specific heat for a collection of 
five needle-like crystals with a total mass m = 4.68(2) 
mg was measured by a long relaxation method utilizing 
a commercial ^He microcalorimeter (Quantum Design 
PPMS-9). In various magnetic fields, the thermometer 
on the calorimeter puck was calibrated before the 
measurements, and the addenda was predetermined in a 
separate run. Magnetic susceptibility measurement was 
performed using a superconducting quantum interference 
device magnetometer (MPMS-5). 


Figure 1 encapsulates the 6 -axis MR (pt) and the 
extracted i?c 2 with fields applied along a*, 6 and c*, 
respectively. To eliminate the ambiguity from the 
superconducting fluctuations, the 50% criterion was used 
in determining i?c 2 , *-e., the field at which pb reaches 
50% of the normal state resistivity. We employed the 
Werthamer-Helfand-Hohenberg (WHH) formula for an 
isotropic one-band BCS superconductor in a dirty limit to 
fit Apparently, at low temperatures, the one- 

band WHH model fails to satisfy the extracted 
By linear extrapolations, i ?“2 j ^c 2 (^) s-nd H ^2 
estimated to be 4.2, 9.6 and 3.3 T, respectively. The 
estimated value of H^ 2 i'^) close to the Pauli-Clogston 
limiting field Hp = 1.84Tc ~ 8.5 T, excluding the 
possibility of triplet pairing. Besides, all of the three 
curves show a linear T-dependence at low temperatures, 
which is consistent with the previous report^^ but in 
sharp contrast with what it should display in most 
superconductorsie., a concave down curvature 
at low temperatures. This anomalous behavior was 
regarded as the evidence for multiband superconductivity 
in iron-based superconductorsMoreover, as shown 
in the inset of Fig. 1(d), the strong T-dependent 
superconducting anisotropy 7 = H^ 2 /^c 2 oc ^c 2 /^c 2 
provides further evidence for multiband scenario as the 
case in two-band superconductor MgB 2 and most iron- 
based superconductors . Using the anisotropic GL 
formula 

17:2(0) = <Ao/ 2 ^C,( 0 )a( 0 ), ( 1 ) 

where (j)Q is the flux quantum and ^j(O) is the GL 
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FIG. 2. (Color online) (a) The low-field magnetization 
curve M{H) under field applied along &-axis at various 
temperatures. The solid red line is fitted to the 1.9 K data, 
(b) The extracted H^iiT) for Ta 4 Pd 3 Tei 6 . The inset shows 
the diamagnetic signal under 10 Oe with zero-field-cooling 
(ZFC) mode. 

coherence length along the j-direction, ^(0) at zero 
temperature is calculated to be 66.1, 151.0 and 51.9 A 
for a*, b and c*, respectively. Obviously, the interchain 
coherence lengths (0) and (0) are much longer 
than the distance between any two adjacent chains in 
Ta 4 Pd 3 Tei 6 denoting the continuous superconducting 
phase across the chains. In other words, the super¬ 
conductivity here is anisotropic but three-dimensional in 
nature, similar to recent discovered QID superconductors 
Nb 2 Pda;Se 5 and Ta 2 Pda,S 5 J ^^i^° For the normal state, the 
longitudinal MR {H || /) response is negligibly small, 
while the in-chain MR (iJ _L i) is obviously large in 
magnitude and positive, which is typical of a QID metal 
due to the change in carrier trajectories induced by the 
Lorentz force 

The lower critical field values, were deter¬ 

mined from low-field magnetization curves M(H) with 
field applied along the 6 -axis as shown in Fig. 2. The 
low-field parts almost overlap with the Meissner fine (the 
solid red line in Fig. 2(a)) due to the Meissner effect. 
Thus, H^i{T) could be defined at the point where M{H) 
deviate by 2% from the perfect Meissner response, and 
the extracted results are plotted in Fig. 2(b). Since 
that our crystal is needle-like and the demagnetization 
factor for field along the needle-like direction is negligibly 
small, H^i{0) can be directly estimated to be 29.9 Oe 
by fitting the extracted data to the formula, Hci{T) = 
i/ci(0)[l — (T/Tc)^], represented by the dashed red line 
in Fig. 2(b). With the results for and 


calculated the GL parameter to be about 84 using 
the equation i?c 2 (0)/i?ci (0) = 2KQL/lnKGL- The above 
results indicate that Ta 4 Pd 3 Tei 6 is an extremely type-II 
superconductor. 

Figure 3 presents the T-dependent specific heat of 
Ta 4 Pd 3 Tei 6 crystals divided by temperature at zero 
magnetic field. A sharp anomaly, denoting the super¬ 
conducting transition, can be clearly observed at Tc ~ 
4.14 K (defined by an entropy conserving construction). 
In Fig. 3(a), the raw heat capacity data between 3.5 
K and 7 K is fitted to the formula C/T = 7 - 1 - /3T^, 
and the fitting gives parameters 7 = 51.2(1) mJ/mol 
and (3 = 13.53(1) mJ/mol K'^, both of which are slightly 
larger than our previous report This small discrepancy 
is possibly due to the measurement precision as well 
as the uncertainty of sample mass weighed. Besides, a 
small residual linear term 70 = 5.1 mJ/mol exists, 
as observed from the intercept of C/T vs in the 
T —^ 0 limit, the origin of which may be the presence of 
nodal quasiparticles for a nodal superconducting gap or 
a small fraction of nonsuperconducting metallic impurity 
phase. By subtracting the phononic contribution, the 
electronic specific heat can be obtained by: Cei(T) = 
C{T) — which is plotted in Fig. 3(b) as Cei/T 
versus T/Tc- The zero-field data at low temperature 
is obviously irreconcilable with the conventional s- 
wave order parameter as the significant quasi-particle 
excitations. 

The so-called a model, which was devised to sim¬ 
ulate the strong-coupling effect, has gained significant 
success in explaining the properties of strong-coupling 
superconductors;^ Within this model, the temperature 
dependence of energy gap A(T) approximately follows 
weak-coupling BCS behavior multiplied by a dimension¬ 
less parameter a. In the BCS theory, the entropy S in 
the superconducting state is given by^ 

o /*27r poo 

S = / / /In/ + [1 - /]ln[l - f]dedp, (2) 

Jo Jo 

where / is the Fermi function / = (1 -|- 
with the quasiparticle energy E = \JT A^(T, </)), 7 „ 
is the normal state 7 for the superconducting part, 
and A(T, (/>) is the temperature and angle dependence 
of the gap function. The electronic heat capacity is 
thus calculated by Cei = T{dS/dT). To elucidate the 
superconducting order parameter, the raw data of Ce\ is 
analyzed and fitted by various models, i.e., the isotropic 
s-wave function, anisotropic d-wave function, two-fold- 
symmetric anisotropic s-wave function [A((/)) = ai + 
0:2 cos 2(j)\ , and two-gap scenarios (s-|-s waves and s-|-d 
waves). The contribution of Ce\ from nodal quasiparticles 
or the small fraction of nonsuperconducting metallic 
impurity phase is regarded as 70 T in our fittings. The 
fitting curves for the cases of single-band d-wave and two- 
band s+d waves are plotted in Fig. 3(b), from which, 
one can unambiguously observe that the two-gap (s-|-d 
waves) model best captures the experimental data from 
the base temperature T/Tc ^ 0.12 up to Tc- The fitting 
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FIG. 3. (Color online) (a) Specific heat of Ta 4 Pd 3 Tei 6 
crystals divided by temperature (C/T) at zero magnetic field. 
The solid pink line is the fit to the normal-state data from 
4.5 to 7 K using the formula C/T = 7 -I- /3T^. (b) T- 

dependent electronic contribution to the specific heat divided 
by temperature (Cei/T). Red and magenta dashed lines are 
the fitting curves with two-gap {s+d waves) and single-gap 
d-wave scenarios. The inset shows the deviation of the fitting 
value for the case of s-l-d waves with the data (represented by 
DF in the plot), (c) expands the low temperature region. 
Lines show the theoretical fitting of various s-wave gap 
functions as described in the text, (d) The deviation of the 
fitting from the data for the cases of s-|-s waves (upper panel) 
and anisotropic s-wave (lower panel). 


gives the following parameters: a® = 0.916 (for s wave), 
= 1.652 (for d wave) and the ratio of weight Jn ■ In = 
33.5 : 66.5, where the partial Sommerfild coefficient 
characterizes each band. Nodal quasiparticles for the 
nodal superconducting gap or impure phase are expected 
to be the source of the finite value 70 in T —>■ 0 limit for 
this case. The fittings with other gap functions (full-gap 
cases) mentioned above deviate the experimental data 
significantly, especially in low-temperature region, the 
close-up view of which is displayed in Fig. 3(c). For the 
full-gap cases, the residual linear term 70 should come 
from a small amount of nonsuperconducting metallic 
impure phase. In order to show more clearly the 
deviations of the fitting curves from the data for the 
case of s-bs waves and anisotropic s-wave, we plots the 
differences in Fig. 3(d). Therefore, our result strongly 
suggests Ta 4 Pd 3 Tei 6 is a two-gap superconductor with 
the pairing symmetry of s+d waves. This fitting 
is consistent with the STS study as the temperature 
dependence of tunneling spectra can be well fitted with 
s+d waves (the anisotropic s-wave symmetry can as well 
describe their data while in our heat-capacity study it 
can not)j^ In addition, by calculating the maximum gap 


FIG. 4. (Golor online) (a) Specific heat of Ta 4 Pd 3 Tei 6 
crystals in selected magnetic fields parallel to c*-axis, plotted 
as C/T vs T^. (b) The plots of C/T vs T^ in the low-T 

region. For each applied field, 7(i7) = limr-jo C{T, H)/T was 
determined from a linear extrapolation of the experimental 
data in the temperature range 0.5—1.3 K. (c) Magnetic field 
dependence of the electronic specific heat coefficient 'y{H) in 
the mixed state. 


value using the derived parameters (A^ax = 0.335*A® + 
o.ees^A®*), we find that A„iax/fcB7c ~ 2.47, very close 
to the value 2.3 reported in Ref. 24. The scenario of 
two-band s-wave gaps with significant ratio has been 
proposed to be possible in understanding the data of 
thermal conductivity!^ However, our analysis of the heat 
capacity does not provide sufficient evidence for this 
possibility. On the contrary, the evidence of anisotropic 
multi-gap structure is consistent with the study of 
thermal conductivity;^ which suggests the presence of 
node in the superconducting gap. Overall, our result 
reported here reconciles the recent experimental studies 
regarding the pairing symmetry in Ta 4 Pd 3 Tei 6 . 

As the magnetic field dependence of specific heat at 
low temperatures is instructive for the quasi-particle 
excitations, we measured the low-T specific heat of 
Ta 4 Pd 3 Tei 6 crystals under magnetic fields {fioH < 
1.4 T) for H II c* as shown in Fig. 4(a). The 
heat capacity anomaly is gradually suppressed to lower 
temperatures with the increase of field. The enlarged 
part of the data is displayed in Fig. 4(b), from 
which, one can observe that minor upturns under 
= 0.3 and 0.5 T appear at low temperatures, 
probably stemming from Schottky anomaly. Therefore, 
we extracted field-dependent electronic coefficient ^{H) 
by a linear extrapolation of 7 (H) = \iu\T^oC{T,H)/T, 
without considering the upturns under 0.3 and 0.5 T. The 
obtained results are presented in Fig. 4(c) as jeiH) = 
7 (iL) — 7 o vs H. In fully gapped conventional type- 
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II superconductors, it is generally known that 7e(-ff) 
should be linearly proportional to the number of field- 
induced vortices,— i.e., ^e{H) oc H. However, we find 
in Ta 4 Pd 3 Tei 6 , a power-law fitting provides the relation 
7 e(i?) j{0A8 jg pig^ 4(c), thereby ruling out 

the isotropic s-wave pairing symmetry. In the case of two 
isotropic s-wave superconducting gaps, ^e{H) is featured 
by an intersection of two straight lines, as observed 
in multiband superconductivity with two s-wave gaps, 
e.g., MgB 2 ,— NbSe 2 ,— and SrPt 2 As 2 .— The nonlinear 
behavior observed here is seemingly in contrast with 
the above scenario. Although ^/~H behavior was once 
considered as the evidence of d-wave superconductors due 
to the Doppler shift of extended quasi-particle excitation 
spectrum,— this kind of nonlinear behavior was also 
frequently observed in many multiband superconductors, 
such as TlNi 2 Se 2 .— Overall, the nonlinear behavior for 
7e(id) could be regarded as an evidence for multiband 
superconductivity. 

In summary, we have investigated the superconduct¬ 
ing properties of Ta 4 Pd 3 Tei 6 crystals. Our measure¬ 
ments indicate Ta 4 Pd 3 Tei 6 is a new anisotropic type- 
II superconductor with the extracted Ginzburg-Landau 


parameter kql = 84. The linear T-dependent Hc 2 at low 
temperatures and the strong T-dependent anisotropy of 
Hc 2 suggest the multiband superconductivity. Indeed, 
the zero-held electronic heat capacity data, Ce\, can be 
satisfactorily described in terms of a two-gap {s+d waves) 
model. The electronic heat capacity coefficient, 7e(i7), 
exhibiting a nonlinear behavior. Thus, Ta 4 Pd 3 Tei 6 
appears to be a rare example of a two-gap superconductor 
with the gap symmetry of s+d waves rather than a multi¬ 
ple fully gapped s-wave symmetry. This superconducting 
state is very interesting and worthy to investigate further, 
especially when noting the QID structural characteristics 
in Ta4Pd3Tei6. 
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